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Thin layers of water on biomolecular and other nanostructured surfaces can be supercooled to 
temperatures not accessible with bulk water. Chen et al. [PNAS 103, 9012 (2006)] suggested that 
anomalies near 220 K observed by quasi-elastic neutron scattering can be explained by a hidden 
critical point of bulk water. Based on more sensitive measurements of water on perdeuterated 
phycocyanin, using the new neutron backscattering spectrometer SPHERES, and an improved data 
analysis, we present results that show no sign of such a fragile-to-strong transition. The inflection 
of the elastic intensity at 220 K has a dynamic origin that is compatible with a calorimetric glass 
transition at 170 K. The temperature dependence of the relaxation times is highly sensitive to data 
evaluation; it can be brought into perfect agreement with the results of other techniques, without 
any anomaly. 

PACS numbers: 87.15.kr,64.70.P-,66.30.jj,61.05.F- 



In contrast to bulk water, protein hydration water can 
be supercooled down to a glass transition at T g ~ 170 K. 
Near T g translational degrees of freedom arrest, which in- 
duces discontinuities in the specific heat and the thermal 
expansion coefficient of the hydration water [1-4]. Due 
to the dynamic nature of the glass transition, freezing of 
microscopic degrees of freedom can already be observed 
far above TL 

The protein dynamic transition is an abrupt onset of 
atomic displacements on the microscopic length and time 
scale probed by quasielastic neutron scattering (QENS). 
First observed twenty years ago in hydrated myoglobin 
and lysozyme at Ta — 240 K |5J], it is now known to 
be a generic property of hydrated proteins, while it is 
absent in dehydrated systems. It is therefore related to 
the dynamics of the hydration shell. 

In QENS, mean squared dispacements (5x 2 ) are de- 
duced from the elastic scattering intensity S(q, 0). Due 
to the finite spectrometer resolution (fwhm 2 Aw), one 
actually measures S(q, \u\ < Aoj). Full spectral mea- 
surements show that the anomalous decrease of the cen- 
tral peak is compensated by increasing inelastic wings 
Q . These effects have been interpreted as precursors 
of the glass transition. Since QENS probes structural 
relaxation at Aw -1 ~ 100 ps, while the calorimetric T g 
refers to a time scale of 100 s, it is natural that Ta is lo- 
cated far above T g : the protein dynamic transition is the 
microscopic manifestation of the glass transition in the 
hydration shell. The time-scale dependence of Ta also 
explains its variation with viscosity [7H9( . 



Recently, these views have been challenged by the 
suggestion that there might be a time-scale indepen- 
dent transition. Support came mainly from QENS ex- 
periments on the backscattering spectrometer HFBS at 
NIST. In hydrated lysozyme, DNA and RNA, a kink was 
found not only in S(q, \u)\ < Aw), but also in the a relax- 
ation time (t) deduced from full spectra S(q,uj) (Tol - |T3l | . 
This change of (r) (T) from high-T super- Arrhenius to 
low-T Arrhenius behavior at Tl ~ 220 K has been inter- 
preted as a fragile-to- strong-transition (FST) from the 
high density (HDL) to the less fluid low density phase 
(LDL) of supercooled water [l4[ . In this view, a qualita- 
tive change in the dynamics occurs when the so-called 
Widom line is crossed, which extrapolates the phase 
boundary beyond the conjectured critical point [la. 1 1 6jj . 

Similar behaviour was found for confined water in var- 
ious non-biological environments (l7l [l8l and references 
therein]; Tl was mostly located between 215 and 228 K. 
There were, however, two remarkable exceptions: mak- 
ing a substrate hydrophobic lowered Tl by 35 K [h|; 
using another neutron spectrometer with wider dynamic 
range lowered Tl by about 10 to 20 K [2(j. Critics ob- 
jected that relaxation times measured by other spectro- 
copic techniques show no anomaly at Tl [2l| - [23l ] . 

In this paper we will argue that the changing temper- 
ature dependence of log((r))is likely to be an artifact of 
the data analysis, and that the kink in the elastic inten- 
sity also admits a conventional, dynamic interpretation. 
We will refer to the hydration shell of different proteins. 
If the anomalies observed around Tl were due to an in- 
trinsic property of water, then the choice of the protein 
should not matter. 
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To improve the empirical base, we measured the mo- 
tion of protein hydration water with enhanced sensitivity, 
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FIG. 1: Neutron scattering spectra of H2O on d-CPC at 
q — 1.0 ± 0.15 A -1 . Note the truncated logarithmic scale. 
Solid lines are fits with (1), convoluted with the instrumen- 
tal resolution (sample data at 100 K). — Inset: central peak, 
normalized to the value at cj = 0. 
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using a new neutron spectrometer and a fully deuterated 
protein. We chose C-phycocyanin (CPC), a light harvest- 
ing, blue copper protein. Deuteration drastically reduces 
scattering from the protein so that we become more sen- 
sitive to weak quasielastic scattering from the hydration 
water. CPC has been investigated by neutron scattering 
before [2^ - |28j . The sample was purified from a prepa- 
ration of Crespi [HI, [2!| with nearly 99 % deuteration. 
The absorbance ratio ^20 nm/^bso nm — 7.3 is above the 
analytic-grade value of 5 [30(. The integrity of the sec- 
ondary structure was ascertained by CD spectroscopy. 
The protein was dialyzed against distilled water, freeze- 
dried and rehydrated to h = 0.3 g/g [27J, the same degree 
of hydatrion as in Ref. 10]. 

QENS was performed on 500 mg of hydrated D-CPC 
with the new backscattering spectrometer SPHERES at 
FRM II [Hj]. It has a resolution (fwhm) of 0.62 fieV 
(HFBS 0.85 /ieV), high flux, and an outstanding signal- 
to-noise ratio of 1000:1. Scattering wavenumbers q range 
from 0.2 to 1.8 A -1 . We chose about the same energy 
range ±10 fieV as customary at HFBS. Spectra at tem- 
peratures T between 100 and 300 K were collected for 6 
to 12 h. An empty cell measurement was subtracted and 
spectra were normalized to a vanadium standard. The 
analysis was performed with the program Frida [32| . 

Fig. Q] displays selected spectra in the relevant temper- 
ature range. The logarithmic intensity scale emphasizes 
both the broad, asymmetric wing of instrumental resolu- 
tion function, unavoidable in a crystal spectrometer, and 
the quasielastic scattering. At least two distinct compo- 



FIG. 2: (a) Scattering intensity integrated over huj < —3 fieV. 
(b) Elastic scattering intensity S(q, 0) and result of (2), com- 
pared to literature data [H, [TTJ] . (c) Thermal expansion co- 
efficient dp of water in hydrated lysozyme powder (+, spec- 
troscopic measurement of hydration water in lysozyme 
solution (x, densitometry, j33j]), of bulk water (line). 



nents of the scattering function S(q,u)) can be observed: 
(i) An elastic component that has the shape of the in- 
strumental resolution function R q {uj). It is mainly due to 
the coherent structure factor of the protein (as evidenced 
by strong small- angle scattering for q <C 1 A -1 as 
well as to a protein incoherent term mainly due to the 
exchanged protons (20%). Also immobilized water near 
charged groups contributes to this component, (ii) A 
quasielastic component that appears above about 170 K. 
It is mainly due to motion of adsorbed H 2 molecules. 

To start with a model-free data characterization, we 
examined the T dependence of the integrated inelastic 
intensity for different u> intervals, preferentially on the 
ijj < side where the resolution wing is weaker. Be- 
tween 180 and 230 K, the inelastic intensity begins to 
rise strongly. There is no sharp kink; the temperature 
of maximum inflection depends on the w range and on q 
(Fig. 2a). This is a first indication for a dynamic phe- 
nomenon rather than a hidden phase transition: the q 
and T dependence is compatible with a smooth onset 
of translational diffusion. The elastic intensity in com- 
pensation shows a step-like decrease at 220 K in paral- 
lel with data obtained with lysozyme. This suggests a 
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common transition temperature of both hydrated pro- 
teins (Fig. 2b). In the following, we concentrate on 
q = 1.0 ±0.15 A" 1 . 

To analyse the spectra, we use the same model as Chen 
et al. [l(| , consisting of an elastic line and a Kohlrausch- 
Williams- Watts function (KWW, Fourier transform of a 
stretched exponential [HI): 



S t h(q, u) = f q \ ai5(u) + a 2 



& iujt-{\t\/rf 



(1) 



The Debye- Waller factor f q = cxp(— (fa 2 ) q 2 ) accounts 
for phonon scattering outside our energy window. The 
mean squared displacement (fa 2 ) is assumed to be linear 
in T; from 100 and 170 K: d{Sx 2 )/dT = 7 • 10~ 4 A 2 /K. 
/? was fixed to 0.5 as in the analysis of Chen et al. 

Equation (1) can fit the data sets for T ^ 240 K, 
hereby covering the temperature range of the proposed 
transition. At higher T, the data start to deviate from 
Eq. (1): additional quasielastic scattering, likely due to 
relaxational or fast anharmonic motion of the protein, is 
not included in this bimodal model. 

In a first attempt to extract the relaxation times, we 
have used a standard fitting procedure for the data eval- 
uation. Fitting the quasielastic wings for T < 240 K, we 
are in the regime ujt 3> 1 where the KWW function goes 
in first order with (wt)~ 1- ' 9 . In such a power law regime, 
the fit parameters amplitude and time scale r are degen- 
erate. Therefore, the 240 K data were used to determine 
the relative contribution of the water 02/(01+02) to 0.42 
which was then used as input for the other temperatures. 

However, this approach neglects a severe limitation of 
the standard fitting procedure for narrow lines: When fit- 
ting the central peak, the convolution of 5th with the ex- 
perimental resolution R q (uj) becomes nontrivial for T < 
220 K where the peak in the KWW function becomes 
so sharp (compared to the experimental energy binning) 
that the straightforward calculation Sth(g, uj) ® R q (uj) as 
a Riemann sum yields a wrong amplitude a\. Because of 
the mentioned degeneracy of amplitude and time scale, 
this results in a distorted (r) . We have therefore repeated 
the fits with an improved fitting procedure, neglecting the 
erroneous a\ and fixing directly 02 at 0.42. 

Fig. 3 shows an Arrhenius plot of the mean relaxation 
time (r) = / dtexp(-t/r)P = /3 _1 r (ft- 1 ) r resulting 
from both fitting procedures. As can be seen, the temper- 
ature dependence of the relaxation times obtained from 
the standard fitting procedure resembles the results of 
Chen et al. This trend is however not observed any more 
in the results of the improved fitting procedure which 
agree well with the structural (a-)relaxation times of wa- 
ter on other proteins taken from literature (also shown). 
There is a good agreement with high-T myoglobin data 
measured by QENS over a wider energy range than cov- 
ered here [| and for T < 220 K with both, 2 NMR [H 
and dielectric spectroscopy Q which show a smooth in- 
crease without inflection at 220 K. The curvature of (r) in 
the Arrhenius plot indicates fragile behaviour according 
to a Vogel-Fulcher relation, (r) = (r) exp [T n /(T - T )}. 
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FIG. 3: Average relaxation times (r) versus reciprocal tem- 
perature for hydration water on different proteins: (■): c- 
phycocyanin with improved fitting procedure and (□): with 
a standard fitting procedure. Error bars were estimated by 
varying P between 0.45 and 0.55. (v): lysozyme, QENS [l(J. 
(A): myoglobin, QENS with wide lu range (xh myo- 
globin, NMR [35fl . (+): myoglobin, dielectric loss [J], ver- 
tical arrow: conjectured Tl of the FST, horizontal arrow: 
instrumental resolution (Aw) -1 . Dashed line: Vogel-Fulcher 
fit to the data above 200 K with fixed T g = 170 K, yielding 
(t) = 6.5 • lO -13 s, T H = 745 K, T = 147 K. 



For low T, our (r) exceeds the resolution Aw -1 ~ 2 ns 
by some orders of magnitude. Still, for \huj\ > 2 /xeV 
there are clear quasielastic wings that allow to determine 
r within the chosen model with fixed parameters 02 and 
P. Therefore, this analysis can be used to probe the 
spectral behaviour for the predicted FST discontinuity. 

Whereas our first analysis, using the standard fit- 
ting procedure, reproduces the results from hydrated 
lysozyme [l(| - the basis for the postulated FST at 220 K 
- we can rule out any discontinuity at this temperature 
using the improved fitting procedure. We suspect there- 
fore that also the various results backing the FST are af- 
fected by the problems of the standard fitting procedure: 
In the power-law regime ujt 3> 1 , fitted (r) are meaning- 
ful only if the amplitude 02 has been fixed at a physically 
reasonable value. In the FST reports, the amplitude, if 
mentioned at all, is just called a scaling constant [17j |. 
The fitted (r) (T) are therefore possibly distorted by an 
uncontrolled drift of 02 (T). This could explain why the 
cross-over temperature Tl varies with pressure, while the 
cross-over time appears constant [l3T | . 

Below approximately 200 K, the relaxation times ex- 
tracted from our improved technique deviate slightly 
from the Vogel-Fulcher relation. Being orders of mag- 
nitude above the instrumental resolution, a detailed dis- 
cussion of this effect seems inadequate. If it were real, 
one could attribute it to the onset of the glass transition 
which is known to be about 30 K broad [1, [36|. The 
observed process would then no longer be the isotropic 
a relaxation 1351. 



4 



In contrast to the kink in (r) , the kink in the observed 
elastic scattering intensity S(q, |w| < Aw) is an undis- 
putable experimental fact. Fig. 2b shows that observa- 
tions at SPHERES (CPC) and HFBS (lysozyme pjj) 
almost coincide below 240 K. The difference at higher T 
can be explained by our use of deuterated CPC, which 
reduces the relative contribution of the protein. 

The expected elastic scattering is 

S cxp {q,0) = J du S th (q,uj) R q {-uj) 

{pOO J j_ \ (2) 

axR g (0) + (l-ai) J — R q (t)> 

with R q (t) — Jdui cos(uit)Rq(u)). Indeed, the low-T data 
are perfectly reproduced by Eq. (2) using the values from 
the Vogel-Fulcher relation for (r) . This confirms the long 
established explanation of the protein dynamic transition 
as a time-scale dependent glass transition: When t _1 
exceeds the instrumental resolution Aw, then the central 
part of the KWW spectrum starts to broaden so that 
S(q, \u\ < Aw) falls below 

This is further corroborated by elastic scattering (q — 
1.4 A -1 ) from hydrated myoglobin on the thermal neu- 
tron backscattering spectrometer IN13 [5j. Its resolu- 
tion width 2 h Aw = 8 fieV is an order of magnitude 



larger than that of HFBS or SPHERES, and the kink 
in S(q, |w| < Aw) is observed accordingly at a higher 
temperature, namely the often reported Ta — 240 K. 

Previous QENS experiments were motivated by the 
hope that hydration layers would provide an opportu- 
nity to study water in a deeply supercooled state that 
is not accessible in the bulk. However, thin hydration 
layers seem to be qualitatively different from bulk water 
[211 ] . This view is also supported by literature data of the 
thermal expansion coefficient ap of lysozyme hydration 
water (Fig. 2c). In concentrated solution, ap is nearly T 
independent and much larger than the bulk value (33j . In 
hydrated powder, ap is estimated from the O-H stretch- 
ing frequency of water, which varies with the average H 
bond length. A striking step is found near 170 K, induced 
by the softening of the O-H-0 hydrogen bond network 
J3|. No anomaly is observed at Tl ~ 220 K. 

All observations agree, however, with a glass transition 
of the hydration shell and a dynamic onset temperature 
Ta, which varies with the probe frequency (Fig. 2b). The 
different values of ap for bulk and hydration water could 
indicate that the protein acts as a patch-breaker [37l ]. 
suppressing critical fluctuations present in bulk water. 

This work has been supported by Deutsche For- 
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